Block Antiferromagnetism and Checkerboard Charge Ordering in Alkah-doped Iron 

Selenides Ri-xFe2-ySe2 



Wei Li/ Shuai Dong,^ Chen Fang,^ and Jiangping Hu'^'^'j^] 

^Department of Physics, Fudan University, Shanghai 200433, China 
^Department of Physics, Southeast University, Nanjing 211189, China 
^Department of Physics, Purdue University, West Lafayette, Indiana ^1901, USA 
^Beijing National Laboratory for Condensed Matter Physics, 
Institute of Physics, Chinese Academy of Sciences, Beijing 100080, China 
(Dated: October 4, 2011) 

By performing first-principles electronic structure calculations and analyzing effective magnetic 
model of alkali-doped iron selenides, we show that the materials without iron vacancies should ap- 
proach a novel checkerboard phase in which each four Fe sites group together in tetragonal structure. 
The checkerboard phase is the ground state with a block antiferromagnetic (AFM) order and a small 
charge density wave order in the absence of superconductivity. Both of them can also coexist with 
superconductivity. The results explain mysterious 2x2 ordered patterns and hidden orders ob- 
served in various different experiments, clarify the missing link between AFM and superconducting 
phases, suggest that the block-AFM state is the parent state, and unify the understanding of various 
observed phases in alkali-doped iron selenides. 

PACS numbers: 74.70.-b, 74.25.Jb, 74.25.Ha, 74.20.Mn 



The newly discovered alkali-doped iron selenide 
superconductors pj-t3] have attracted much research at- 
tention because of several distinct characters that are 
noticeably absent in other iron-based superconductors, 
such as the absence of hole pockets at F point of Bril- 
louin zone in their superconducting (SC) phases [4^,6^ and 
AFM ordered insulating phases [7]-[9] with very high Neel 
transition temperatures in their parental compounds [3 . 
Due to these distinct physical characters from their pnic- 
tide counterparts, Ri-x^e2-ySe2 are expected to be 
ideal grounds to test theoretical models of iron-based 
superconductors. Models based on different mecha- 
nisms have suggested different pairing symmetries for 
Ri-xFe2-ySe2: weak coupling approaches based on 
spin-excitation mediated pairing predict a d-wave pair- 
ing symmetry [TOH13], strong coupling approaches p!4] - [T6] 
which emphasize the importance of next nearest neighbor 
(NNN) AFM local exchange coupling suggest the pair- 
ing symmetry is a robust 5- wave, not different from the 
5'^-wave symmetry obtained in their pnictide counter- 
parts and models with orbital fluctuation mediated pair- 
ing suggest a ^^^-wave pairing for both iron selenide and 
pnictide materials [17]. 

While the iron-selenide superconductors have gener- 
ated considerable excitement, there are deep confusions 
regarding the delicate interplay between Fe vacancies, 
magnetism and superconductivity. Many latest experi- 
mental results in Ri_x^e2-ySe2 indicate that the insu- 
lating AFM and SC phases are phase separated[18-23 . 
In particular, the recent scanning tunneling microscopy 
(STM) measurements on Ki_a;Fe2-2ySe2 clearly suggest 
phase separation [21] 122 • The material was shown to 
be phase separated into iron vacancy ordered regions 
and iron vacancy free regions. The former is insu- 



lating and shows a a/5 x a/5 vacancy ordered pattern 
while the latter is SC. The neutron scattering experi- 
ments have shown that the ^/b x a/5 vacancy ordered 
phase is also AFM ordered [7 . More surprisingly, the 
AFM order is affected by SC pairing [8 , a result diffi- 
cult to be understood within the picture of phase sep- 
aration. Besides the a/5 x a/5 vacancy ordered phase, 
an additional insulating phase with a 2 x 2 ordered pat- 
tern was also observed [19| |20l [24]. Moreover, in the 
SC state where there is few vacancy, both STM [211 122] 
and angle resolved photoemission spectroscopy (ARPES) 
experiments [6 suggest that there is an additional sym- 
metry breaking order in the SC phases Ell ED- The 
microscopic origin of this order and how it is related to 
the AFM phase are not understood. 

In this Letter, we show that alkali-doped iron selenide 
superconductors without iron vacancies should approach 
a checkerboard phase in which each four Fe sites group 
together in a tetragonal structure. This broken symme- 
try state is essentially driven by the same magnetic ex- 
change couplings that drive the insulating AFM phase in 
the a/5 X a/5 vacancy ordered state. We perform first- 
principles electronic structure calculations and develop 
an effective magnetic model to show the existence of such 
a broken symmetry state. The checkerboard phase is the 
ground state with a block-AFM (BAF) order in the ab- 
sence of superconductivity. The BAF fluctuations and 
the checkerboard lattice distortion are strongly coupled. 
A weak BAF order and the checkerboard lattice distor- 
tion can coexist with superconductivity. These results 
essentially suggest the BAF ordered state are the par- 
ent state of alkali-doped iron selenide superconductors. 
The results consistently explain the mysterious 2 x 2 or- 
dered pattern which was misunderstood as another va- 
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FIG. 1: (Color online) (a) The crystal structure of KFe2Se2 
compound in our calculation. The system consists of K 
(royal), Fe (dark yellow) and Se (green) atoms; (b) Schematic 
checkerboard lattice structure and spin ordering pattern in 
the BAF state. The lattice distortion is labeled by lattice 
constants X and X' . The magnetic exchange couplings are 
also indicated. 



cancy ordered state and the STM and ARPES exper- 
imental results [6, 21 . This study clarifies the missing 
link between AFM and SC phases and essentially unifies 
the understanding of various observed phases [IP. 

Results from LDA calculations We start with the fol- 
lowing simple question: if the system is free of iron 
vacancies, what should be the ground state if it does 
not become SC? To answer this question, we perform 
the first principles calculation to investigate the ground 
state of an iron vacancy free domain. We calculated 
the energy of a number of different possible magneti- 
cally ordered states, including non-magnetic (NM), ferro- 
magnetic (FM), collinear-AFM (CAP, the state observed 
in iron-pnictides[25^ 24), bicollinear-APM (BCAP, the 
state observed in PeTe) [27H3Q and BAP whose pattern is 
shown in Pig. [ijb) where four Pe sites group to form a su- 
per cell. All these calculations were performed using the 
projected augmented wave method [32] as implemented 
in the VASP code [33 , and the Perdew-Burke-Ernzerhof 
(PBE) exchange correlation potential [34 was used. A 
500eV cutoff in the plane wave expansion ensures the 
calculations converge to 10"^ eV. For the BAP state, all 
atomic positions and the lattice constants were optimized 
until the largest force on any atom was 0.005eV/A. We 
used a 9 X 9 X 5 Monkhorst-Pack k-grid Brillouin zone 
sampling throughout all of calculations. 

The ground state energies, the magnetically ordered 
moment and their lattice constants in various states are 
listed in Table [ij The BAP state clearly has the low- 
est energy. Without turning on /7, our LDA results 
show that the BAF ordered state is metallic. However 
if a finite U > l.OeV is added, the LDA+/7 calcula- 
tion shows that the state becomes insulating. There are 
strong lattice distortion in the BAF state. The lattice 
constant X (between two nearest sites in one supercell) 



TABLE L Geometric, energetic and magnetic properties of 
KFesSes. Results in the NM/FM/AFM/CAF/BCAF/BAF 
and optimized BAF configurations using fully optimized 
structures are all shown. is the total energy difference 

per iron atom in reference to the unoptimized experimental 
structure[l], and mpe is the local magnetic moment on Fe. 



KFe2Se2 


AE (eV/Pe) 


a(A) 


c(A) 


mFe(Ms) 


NM 





3.9136 


14.0367 





FM 


-0.2400 


3.9136 


14.0367 


2.781 


AFM 


-0.2384 


3.9136 


14.0367 


2.135 


CAP 


-0.3510 


3.9136 


14.0367 


2.446 


BCAP 


-0.3159 


3.9136 


14.0367 


2.556 


BAF 


-0.3127 


3.9136 


14.0367 


2.552 


BAF(opt.) 


-0.3568 


3.8553 


14.4099 


2.635 



and X' (between two supercells) as labeled in Pig. [TJb) 
are X=2.6A, X'=2.85A for U = OeV and X=2.59A, 
X'=3.03A for U = 2eV. This lattice distortion is com- 
parable to the lattice distortion in the ^/5 x ^/5 vacancy 
ordered K0.8Pe1.6Se2 phase[7]. 

Prom Table [l[ the energies of the BCAP state and the 
BAP state are almost degenerate in the unoptimized lat- 
tice tetragonal structure. The actual ground state, in 
fact, depends on the optimization of lattice structure. 
In the previous first principle calculations, the BCAP 
state was shown to be the lowest energy state in the op- 
timized lattice which has monoclinic distortion [35] l36] . 
The monoclinic distortion and the BCAP order can be 
strongly coupled with each other because they break the 
same lattice rotational symmetry. Such a strong cou- 
pling was observed in PeTe in which a single strong first 
order phase transition where both the BCAP ordering 
and the tetragonal- monoclinic distortion take place [27]. 
However, in KPe2Se2, no monoclinic lattice distortion has 
been observed. Without allowing the monoclinic lattice 
distortion, the BAP state becomes the ground state. 

Prom the similar symmetry analysis, the BAP state 
must strongly couple to the lattice distortion shown in 
Pig. [ijb). The lattice distortion quadroples the lat- 
tice unit cell to form a checkerboard pattern. In such a 
checkerboard lattice, the charge ordering can take place. 
If we calculate the electron density distribution around 
Se atoms in the Se-layer above the iron-layer, a charge 
ordering on the Se-layer is observed as shown in Pig. [2j 
which was observed in recent STM experiment [2T] . 

Magnetic models Now we discuss the effective magnetic 
model that can interpret above calculation results. It has 
been shown a magnetic exchange J1-J2-J3-K model [38] 
where Ji, J2 and J3 are the nearest neighbor (NN), NNN, 
and the next NNN (NNNN) exchange couplings respec- 
tively, and is a spin biquadratic coupling term between 
two nearest neighbor sites is a good approximation to 
describe iron-chalcogenides when the lattice distortion is 



FIG. 2: (Color online) Charge density distribution (e/bohr^) 
in the (001) plane crossing the first Se atoms layer with BAF 
order state within the LDA+t/ (U — 2eV) calculations. 



ignored f29l [30l |37l [38]. In Table |T) the energies of the 
BCAF and BAF states in the unoptimized lattice tetrag- 
onal structure are almost degenerate. This degeneracy 
is a strong support of the model since these two states 
are exactly degenerate in the J1-J2-JS-K classical spin 
model[38]. In the magnetically ordered state, the lat- 
tice distortion takes place and the tetragonal symmetry 
is broken. The nearest neighbor exchange coupling Ji 
can take two different values Ji and J( as shown in Fig. 
[TJb). The biquadratic coupling K can be decoupled and 
treated as an effective difference between Ji and J[ as 
well[38l [39]. In general, the NNN J2 can also take two 
different values, J2 and J2 as also shown in Fig. [ijb). 
However, as being proved in other iron-based supercon- 
ductors, the NNN coupling J2 is rather robust against 
lattice distortion. The difference between J2 and J2 is 
rather small. Therefore, the effective magnetic exchange 
model in magnetically ordered state is given by Ji- J{- J2- 
J3 with Ji being strongly FM and J2,3 both being AFM. 
J[ can be weak FM or weak AFM. The similar model has 
been shown to describe the magnetism of the ^/5 x 
vacancy ordered K0.8Fe1.6Se2 phase[9l El]- Therefore, 
while the exact values of the magnetic exchange couplings 
can not be accurately obtained from LDA calculations, 
since the lattice distortions in both cases are similar, it 
is reasonable to believe that these values should not be 
too different from those of K0.8Fe1.6Se2 which has been 
measured by fitting neutron scattering experiments [9 . 
The measured values, which are specified in Fig. |3) 
also give the BAF order ground state. The saved en- 
ergy from magnetic exchange coupling per site is given 
by (— Ji + 2J3 + J[)S'^ . This energy is sightly smaller 
than the saved magnetic energy in the vacancy ordered 
Ko.8Fei.6Se2^ [37]. The spin wave dispersion and the 
imaginary part of dynamic spin susceptibility of the BAF 
phase are shown in Fig. [3| 

There is another interesting prediction if the same 
magnetic model describes both K0.8Fe1.6Se2 and 
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FIG. 3: (Color online) Spin wave dispersion (up) and imagi- 
nary part of dynamic susceptibility (down) for the BAF state 
with Ji-J[-J2-J3 model at Ji = — 36meV, J[ = 15meV, 
J2 = 14meV, J3 = 9meV, taken from Ref. 9 The profile 
of the imaginary part of the dynamic susceptibility is plotted 
at various energies with an energy resolution of 5meV, and it 
is given in arbitrary unit. 



Ka^Fe2Se2. As mentioned before, K0.8Fe1.6Se2 and 
Ka;Fe2Se2 are two phase separated regions. If the same 
effective magnetic model describes both structures, it is 
very interesting to inquire into the magnetic configura- 
tions near the boundary of these two structures. We per- 
form a Monte Carlo (MC) simulation on the J1-J2-J3-K 
model [38 to address this problem. A simple numerical 
simulation, which includes a standard Markov Chain MC 
simulation followed by a zero-temperature relaxation pro- 
cess, is performed to qualitatively investigate the mag- 
netic orders near phase boundaries. A two-dimensional 
spin lattice [Lx x {Lyi-\-Ly2)] is used with periodic bound- 
ary conditions (PBCs). Vacancies with the ^/5 x 
pattern is created in Lyi regions for the K0.8Fe1.6Se2 
phase. A general result we obtained as shown in Fig. 
|4] is that the spin directions between Lyi and Ly2 re- 
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(a) MC snapshot 
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(b) after relaxation 
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FIG. 4: (Color online) (a) A typical MC snapshot (after 3 x 
10^ MC steps) of the classical J1-J2-J3-K magnetic model, 
(b) The spin pattern after the zero-T relaxation. In (a) and 
(b), spins in regions without vacancies and with vacancies 
are in blue and red, respectively. Black circles denote the 
vacancies. 



gions are noncollinear. Since experimentally, the ordered 
AFM moment is along c-axis in K0.8Fe1.6Se2 [7], this re- 
sult suggests that the ordered moment in the BAF state 
must be in the plane. This non-collinearity stems from 
the presence of vacancies and intrinsic magnetic frustra- 
tion among the magnetic exchange couplings, similar to 
the study in the frustrated J1-J2 model [40 . Recent STM 
results have provided an evidence supporting this predic- 
tion. It was shown that the magnetic moment induced 
by individual vacancy in the SC state is indeed in the 
plane [22|. 

The phase separation between the vacancy ordered 
BAF state and the SC state has blurred the interplay be- 
tween magnetism and superconductivity in alkali-doped 
iron selenide. The above results also clarify the connec- 
tion. In iron-pnictides, as increasing doping suppresses 
CAF order, superconductivity develops. The magnetic 
order is able to coexist with superconductivity in a small 
doping region [41 . Even in the region where the magnetic 
order is completely suppressed, orthorhombic lattice dis- 
tortion which couples the fluctuating short range CAF 
order[42j 03], can survive and coexist with superconduc- 
tivity. Our result suggests that the similar physics can 
take place in i?i_a^Fe2-ySe2. The absence of iron vacan- 
cies in the SC state suggests that the true SC material 
has a chemical formula i?i_a;Fe2Se2. The parent state of 
this material should be a BAF state. Increasing doping 
suppresses the BAF state and leads to SC. While it is still 
difficult to determine whether the BAF and SC can co- 
exist, we can safely argue that, similar to iron-pnictides, 
a lattice distortion as shown in Fig. [ijb) that couples to 
the short range BAF fluctuation should be able to coexist 
with SC. 

This picture provide explanations to many puz- 



zling phenomena observed in alkali-doped iron selenides 
i?i_a;Fe2-2ySe2. First, in ARPES measurements, a weak 
but large electron pocket at F point was observed [6]. This 
pocket is almost identical to the electron pockets at M 
point, suggesting the electron pocket is a folded pocket 
due to translational symmetry breaking in the SC state. 
Moreover, in recent STM experiment , a 2 x 2 charge 
density modulation with respect to Fe-lattice was ob- 
served to coexist with the SC phase. These electronic su- 
perstructure are consistent with the checkerboard phase. 
Second, neutron scattering experiments suggested that 
the vacancy ordered-AFM state interacts strongly with 
superconductivity [8 . In a phase separation scenario, 
such a strong interaction is hard to understood. Our 
results resolve such a dilemma. The experiment can be 
easily understood because the BAF state strongly inter- 
acts with both the vacancy ordered-AFM state and the 
SC state. The development of superconductivity is ex- 
pected to strongly suppress the BAF state. Finally, the 
2x2 ordering in insulating samples observed by transmis- 
sion electron microscopy (TEM) [19l [24] can be naturally 
interpreted as our checkerboard state with the BAF or- 
der. Previously it was interpreted as another vacancy 
ordered phase p!9l [24]. Such an interpretation is unlikely 
because the TEM signal of this order is much weaker 
than that of the a/5 x ^/S vacancy order. 

In summary, we show that alkali-doped iron selenides 
Ri-x^e2-ySe2 has a checkerboard phase in which each 
four Fe sites group together in a tetragonal structure. 
The checkerboard phase approaches a BAF order in the 
absence of superconductivity. The phase also exhibits 
small charge density modulation on Se sites. Magnetic 
properties related to this state are calculated. Combining 
with the strong experimental evidence of phase separa- 
tion between vacancy ordered and vacancy free phases, 
we suggest the checkerboard phase is the parent state of 
the superconductor. 
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